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Abstract. We present sensitive, high-resolution radio ob- 
servations of the circumnuclear region of the barred spiral 
galaxy NGC6951. These observations reveal a ring of radio 
emission with many discrete components and a marginally 
resolved nuclear component. We compare the radio ring 
with observations at other wavelengths, and discuss the 
nature of the compact radio components. 
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1. Introduction 

The dynamics and relationships between gas inflow pro- 
cesses, circumnuclear star formation and the presence of 
an active nucleus are important, since they might pro- 
vide valuable leads towards understanding the transport 
of gaseous material towards the nuclear regions of ac- 
tive galaxies. An interesting class of objects for study- 
ing these aspects are the Sersic-Pastoriza or S-P galaxies 
(Sersic & Pastoriza 1965; Sersic 1973) which include galax- 
ies with diffuse and amorphous nuclei in addition to the 
well-known hot-spot systems. Sersic & Pastoriza (1967) 
found that galaxies with these unusual nuclear morpholo- 
gies occurred in barred or mixed-type galaxies. The S-P 
galaxies reflect a broad spectrum of properties with some 
of them also harbouring an active nucleus. The existence 
of both starburst and an active nucleus makes them par- 
ticularly interesting for studying their kinematics, evolu- 
tion and possible relationships between the two forms of 
activity. Circumnuclear star formation has been seen in 
many hot-spot galaxies, and there have been suggestions 
of a correlation between circumnuclear rings and nuclear 
activity (cf. Arsenault 1989). 

Theoretical studies suggest that circumnuclear rings 
arise due to a bar-driven inflow of gas and dust to an in- 



ner Lindblad resonance (ILR) or between two ILRs (e.g. 
Combes & Germ 1985; Athanassoula 1992; Byrd et al. 
1994; Piner et al. 1995). The dense gas which accumu- 
lates in an ILR ring leads to a high star-formation rate 
either due to collisions of the molecular clouds (Combes 
& Gerin 1985), or gravitational collapse in the ring when 
the density reaches a critical value (Elmegreen 1994). To 
understand the formation and fuelling of the AGN, the 
gas must flow inwards from the ILR to sub-parsec scales, 
and this is at present not well understood (cf. Axon & 
Robinson 1996). Some simulations suggest small, steady 
inflow (Piner et al. 1995), while others suggest rapid inflow 
for brief periods (Wada & Habe 1992; Heller & Shlosman 
1994). 

In this paper we report high-resolution radio observa- 
tions of the nuclear region in the barred, late-type grand- 
design galaxy NGC6951, which we had earlier observed as 
part of a survey of 47 S-P galaxies with the Very Large 
Array (VLA) at A20 and 6 cm (Saikia et al. 1994). The 
observations reported here are of higher resolution and 
better sensitivity, and clarify the radio structure of the 
circumnuclear region of this well-studied galaxy. New su- 
pernovae have been reported in this galaxy during 1999 
and 2000 (Cao et al. 1999; Valentini et al. 2000). 



2. NGC6951 

NGC6951, which has been classified as a spiral galaxy of 
type SAB(rs)bc (de Vaucouleurs et al. 1991), has an active 
nucleus (Boer & Schulz 1993; Ho et al. 1995) and a cir- 
cumnuclar ring of star formation (Buta & Crocker 1993; 
Barth et al. 1995). The active nucleus has been classi- 
fied either as a Seyfert 2 (e.g. Boer & Schulz 1993; Ho et 
al. 1995, 1997a) or a LINER (Filippenko & Sargent 1985; 
Munoz-Tuhon et al. 1989; Marquez & Moles 1993). Perez 
et al. (2000) have considered the available data and have 
classified it as a transition object between a 'very high 



D.J. Saikia et al.: The circumnuclear ring in NGC6951 



NGC6951 



8439.9 MHz 



NGC6951 



8439.9 MHz 




I I I I I L 

D 36 38.0 37.5 37.0 36.5 36.0 

RIGHT ASCENSION (B1950) 
Cont peak flux = 9.4073E-04 Jy/beam 
Levs = 5.000E-05 * (-1, 1, 1.500, 2, 2.500, 3, 4, 
5, 6, 7, 8, 9,10,12,14,16,18, 20) 




a 



O 



20 36 37.5 37.0 36.5 36.0 

RIGHT ASCENSION (B1950) 
Com peak flux = 5.3600E-04 Jy/beam 
Levs = 5.000E-05 * (-1 , 1 , 1 .500, 2, 2.500, 3, 4, 

5, 6, 7, 8, 9, 10) 



Fig. 1. Radio images of the nuclear region of NGC6951 at 8.4 GHz with an angular resolution of 2.57 x 1.44 arcsec 2 
along PA 0° (left panel), and 0.86 x 0.71 arcsec 2 along PA 113° (right panel). The peak brightness and contour levels 
are indicated below each image. 



excitation LINER and a possible nitrogen over-abundant 
Seyfert 2.' 

It is at a distance, d, of 24.1 Mpc (Tully 1988), so that 
1" corresponds to 117 pc. Marquez & Moles (1993) have 
noted its high degree of isolation with no other galaxy 
within a projected distance of 1 Mpc and 500 km s _1 in 
redshift, suggesting that it has been free of any gravita- 
tional influence for about 10 9 yr. 

There have been a large number of studies of this 
galaxy which include broad- and narrow-band observa- 
tions at optical and infrared wavelengths (Buta & Crocker 
1993; Marquez & Moles 1993; Barth et al. 1995; Woz- 
niak et al. 1995; Elmegreen et al. 1996; Friedli et al. 1996; 
Knapen 1996; Rozas et al. 1996a, b; Gonzalez Delgado et 
al. 1997; Gonzalez Delgado & Perez 1997; Mulchaey et 
al. 1997; Elmegreen et al. 1999; Perez et al. 2000), spec- 
troscopy of the nuclear region (Filippcnko & Sargent 1985; 
Muhoz-Tuhon et al. 1989; Boer & Schulz 1993; Marquez 
& Moles 1993; Ho et al. 1995, 1997a,b; Perez et al. 2000), 
radio images of the nuclear region (Vila et al. 1990; Saikia 
et al. 1994; Ho & Ulvestad 2001, hereinafter referred to 
as HU) and high-resolution images of the molecular gas 
(Kenney et al. 1992; Kohno et al. 1999). 

The most prominent features are the bar, the bulge 
and the spiral arms, with the bar being better defined 
at infrared wavelengths. The circumnuclear region has a 
pseudo outer, pseudo inner and circumnuclear ring, with 
evidence of isophotal twisting between the bar and the 
circumnuclear ring. The dust lanes appear as straight fea- 
tures from the bar edges. The ring morphology has many 



bright knots, with the largest concentration being towards 
the north-western and south-eastern sides of the ring. The 
hotspots or knots of emission are also seen in the infrared 
K-band and there is a general correspondence between the 
positions of the hotspots in Ha and in K band, although in 
some cases there is some offset in the positions (Knapen 
1996). The CO (1—0) emission shows a twin-peak mor- 
phology which could be due to orbit crowding near the 
outer ILR (Kenney et al. 1992; Kohno et al. 1999). The 
HCN observations by Kohno et al. show this dense gas 
to also have twin-peaks, but these are offset from the CO 
peaks and appear to be closely related to the radio peaks 
in the lower-resolution images. 

3. The nuclear radio ring in NGC6951 

In Fig. 1 we present the convolved and full-resolution radio 
images of NGC6951 at 8.4 GHz. These observations were 
made with the VLA B-array on 1995 October 30 and have 
an angular resolution of 2.57 x 1.44 arcsec 2 along position 
angle (PA) 0°, and 0.86 x 0.71 arcsec 2 along PA 113° 
respectively. The rms noise level in each of the images is 
^0.02 mJy/beam. 

The low-resolution image shows two prominent outer 
peaks at RA 20^ 36 m 36. s 10, Declination 65° 55' 49."00, 
and RA 20' 1 36 m 37/20, Declination 65° 55' 44."25, with 
peak flux densities of 0.94 and 0.86 mJy/beam respec- 
tively. The peak flux density of the central component 
is 0.81 mJy/beam. The high-resolution image shows the 
beautiful ring oriented at a position angle of ^135°, with 
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Table 1. The compact components 



Source 




Sgauss 


Lum 

xlO 19 




mjy/b 


mjy/b 


W/Hz 


36.02+49.1 


0.234 


0.23 


1.60 


36.05+46.3 


0.116 


0.10 


0.69 


36.07+47.8 


0.205 


0.20 


1.39 


36.16+45.2 


0.125 


0.07 


0.49 


36.19+50.0 


0.270 


0.24 


1.67 


36.44+43.6 


0.097 


0.08 


0.56 


36.44+50.6 


0.157 


0.14 


0.97 


36.64+46.3™ 


0.536 


0.54 


3.75 


36.82+49.2 


0.105 


0.09 


0.63 


36.96+42.5 


0.109 




(0.76) 


36.98+48.4 


0.128 


0.11 


0.76 


37.08+42.3 


0.145 




(1.01) 


37.12+45.8 


0.166 




(1.15) 


37.16+43.0 


0.202 




(1.40) 


37.20+44.7 


0.335 


0.30 


2.08 



individual peaks of radio emission. The extent of the radio 
ring is ^9.1x6.2 arcsec 2 , which corresponds to a linear size 
of 1060x720 pc 2 . The positions and flux densities of the 
different components were estimated by fitting Gaussians 
and an underlying tipped plane at the peaks of radio emis- 
sion which are at least ^5 times higher than the rms noise 
level. These values are subject to significant uncertainties 
because of difficulties in defining the background flux den- 
sity levels in the vicinity of each of the knots. Due to the 
irregular background and limited angular resolution of the 
observations we could not get satisfactory fits for four of 
the knots. All the components are listed in Table 1, along 
with their flux density estimates in units of mJy/beam. 
The components have been named such that the first four 
digits represent the seconds of right ascension and the dig- 
its after the + sign the arsec of declination. We have listed 
the peak flux densities of the components from the image 
(column 2), as well as those obtained from the Gaussian 
fits (column 3). This gives an idea of the contribution of 
the background at the different locations. The radio lumi- 
nosities of the individual components are listed in units of 
10 19 W Hz -1 using the flux densities from the Gaussian 
fits, if available. The luminosities of those without satis- 
factory Gaussian fits have been estimated from the peak 
values in the image and are enclosed in parentheses. The 
nuclear component has been marked with the superscript 
n. 

The size and orientation of the ring seen at radio fre- 
quencies is roughly similar to that seen at optical, Ha, 
infrared and mm wavelengths. The correspondence be- 
tween the positions of the radio peaks and those seen at 
other wavelengths is not precise although they occur in 
the same general region. For example, by aligning the nu- 
clear components, it can be seen that in the north-western 
region, the radio knot 36.44+50.6 appears spatially coin- 
cident with a bright peak in the B 1 — K 1 image (Perez et 



al. 2000) as well as a prominent peak in the Ha contours 
(Gonzalez Delgado et al. 1997). On the other hand, there 
are no such corresponding bright features for the nearby 
brighter radio components 36.19+50.0 and 36.02+49.1. 

The separation between the successive peaks of radio 
emission listed in Table 1 ranges from ~0."8 to 3. "5. These 
separations are reasonably uniformly distributed between 
0."8 and 2. "6, with only the largest separation being be- 
yond this range. There does not appear to be a regular 
period in the separation of the peaks of emission. If the ra- 
dio peaks trace the star forming regions, one would expect 
these to be regularly spaced in the gravitational instability 
models (Elmegreen 1994). 

The two most prominent regions of radio emission are 
towards the north-west and the south-east. From CO (1-0) 
and HCN (1-0) observations, Kohno et al. (1999) find that 
although both CO and HCN emission are dominated by 
'twin-peaks' morphology, the HCN emission is better cor- 
related spatially with the peaks of radio emission, which 
are identified with the regions of massive star formation. 
In this case, the peaks of CO emission are identified with 
the region of crowding of the x\ and x-i orbits. 

3.1. Nature of radio emission 

The radio emission from the circumnuclear region is pos- 
sibly a mixture of thermal and non-thermal emission with 
the discrete components being either supernova remnants 
(SNRs) or Hn regions. The diagnostics which are usually 
employed to distinguish them are the brightness temper- 
atures, as in the starburst galaxy Arp 220 (Smith et al. 
1998), and radio spectral indices. The radio components 
are usually too weak for polarization measurements. 

The brightness temperature of the components in 
NGC6951 is small, < 15° K, and does not provide an un- 
ambiguous resolution of the nature of the compact compo- 
nents. Spectral index information on the individual knots 
is not available; however the spectral index, a defined as 
Soc v^ a 1 of the entire circumnuclear emission estimated 
from scaled-array VLA observations between A20 and 6 
cm is 0.84 (HU). Convolving the A3. 6 cm image to lower 
resolutions and comparing with the images of Saikia et al. 
(1994) and HU also suggests a non-thermal spectrum. 

However, a few caveats need to be borne in mind while 
examining the nature of the emission of the components 
based on their spectra. With inadequate resolution H n 
regions and SNRs may not be resolved from each other 
(cf. Gordon et al. 1999). Also, although the overall spec- 
trum from coarser resolution might appear non-thermal, 
individual compact components identified from higher- 
resolution observations may be due to H n regions. An ex- 
ample of this is the blue compact galaxy Henize 2-10 which 
is an archetypal starburst galaxy containing a large pop- 
ulation of Wolf-Rayet stars (Conti 1991). The integrated 
radio spectrum estimated from observations largely with 
an angular resolution of about a kpc has a non-thermal 
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spectrum, with a ~ 0.54, while the compact radio knots 
observed with a resolution of ^30 pc have a flat or inverted 
spectrum (Kobulnicky & Johnson 1999). These are iden- 
tified with dense, compact H n regions which are strong 
sources of free-free emission. For example, the features in 
Henize 2-10 (d~9 Mpc) have a luminosity of ~9xl0 18 W 
Hz -1 at A3. 6 cm, electron densities of up to 5000 cm~ 3 and 
harbour the youngest 'super-star clusters' (Kobulnicky & 
Johnson 1999). A similar feature has been reported in 
the nuclear region of NGC5253 (Turner ct al. 1998, 2000) 
which has a luminosity of ^2xl0 19 W Hz -1 at A2 cm. In- 
frared observations confirm it to be an H II region excited 
by a dense cluster of young stars (Gorjian ct al. 2001). 

The nature of the compact radio components in other 
galaxies with either a circumnuclear ring or a nuclear 
starburst might provide some clues towards understand- 
ing the nature of the components in NGC6951. A num- 
ber of these galaxies have been observed at more than 
one frequency. Amongst the galaxies with radio rings, the 
radio spectrum of the ring in NGC613 was found to be 
~ 0.65 ±0.10 (Hummel & Jorsater 1992), the components 
in NGC1097 (Hummel et al. 1987) and NGC1365 (Saikia 
et al. 1994; Sandqvist et al. 1995; Stevens et al. 1999) 
tend to have non-thermal spectra, as well as the ring of 
emission in the Seyfert galaxy NGC7469 (Wilson ct al. 
1991). In NGC4736, which has a continuum ring, Duric & 
Dittmar (1988) found 9 of their compact components to 
have a thermal spectrum between A20 and 6 cm, and the 
remaining 10 to have a non-thermal spectrum. Turner & 
Ho (1994) observed this galaxy at A6 and 2 cm, and found 
9 of their 14 sources to have a < 0.2. 

The compact components in these systems are a mix- 
ture of both SNRs and H II regions occuring in an area of 
high star formation. The proportion of these two types of 
objects has been investigated in some detail in the archety- 
pal starburst galaxies M82 (Unger et al. 1984; Kronberg 
et al. 1985) and NGC253 (Antonucci & Ulvestad 1988). 
Spectral studies of selected components in M82 (Wills ct 
al. 1997; Allen & Kronberg 1998) showed that of the 26 
components, only one was identified with a thermal emis- 
sion spectrum. However, more recent observations of a fur- 
ther 20 largely weaker components suggest that only 5 of 
these are SNRs while the remaining 15 are likely to be 
Hn regions (McDonald 2001). Thus about two-thirds of 
the compact objects are identified with SNRs. In NGC 
253 about a half of the sources with spectral informa- 
tion have flat spectra and are dominated by emission from 
Hn regions, while the remaining half have steep spectra 
and are identified with SNRs (Ulvestad & Antonucci 1991, 
1997). A similar situation is also seen in the nearby merger 
NGC4038/4039 (the Antennae), where about a third of 
the objects have been suggested to be thermal and the re- 
maining two-thirds non-thermal (Neff & Ulvestad 2000). 
The spectral indices of the components in NGC1808 tend 
to be nonthermal (Saikia et al. 1990; Collison et al. 1994). 
In the starburst galaxy NGC2146 which has an S-shaped 



structure with several knots of emission (Kronberg & Bier- 
mann 1981; Saikia et al. 1994), spectral index measure- 
ments suggest three sources to be good candidates for 
SNRs or radio supernovae, and 6 to be associated with 
compact, dense Hn regions (Tarchi et al. 2000). Based 
on these studies of a wide variety of starburst galaxies, it 
is reasonable to assume that the ring in NGC6951 con- 
sists of a mixture of both these types of objects, although 
our linear resolution of about 100x83 pc 2 is inadequate 
to identify individual SNRs and H n regions. 

3.2. Radio luminosity and supernova rate 

The compact components seen in NGC6951 have a me- 
dian luminosity of 10 19 W Hz -1 at A3. 6 cm. The value 
is similar when one considers those with Gaussian fits as 
well as for the entire list. The luminosity distribution is 
skewed towards lower values, with a tail extending towards 
higher luminosities, suggesting that we have not been able 
to identify the weaker components from the existing ob- 
servations. 

It is interesting to compare the values obtained for 
NGC6951 with those obtained for other galaxies exhibit- 
ing high star formation rates. For ease of comparison we 
have considered those for which observations have been 
made at A3. 6 cm. The median luminosity of the com- 
ponents in NGC253 (d^2.5 Mpc) observed with a linear 
resolution of —3 pc is ~ 4.5xl0 17 W Hz -1 (Ulvestad & 
Antonucci 1991), while those in M82 (d-3.4 Mpc) ob- 
served with a similar linear resolution has a median value 
of ~2x 10 18 W Hz" 1 (Huang et al. 1994; Allen & Kronberg 
1998). Considering more distant galaxies, the median lu- 
minosity for the archetypal starburst galaxy in the south- 
ern hemisphere, NGC1808 (d-10.9 Mpc), is about 8xl0 18 
W Hz -1 with a linear resolution of 30 pc (Collison et al. 
1994), while the knots in NGC1365 (d~20 Mpc) has a me- 
dian luminosity of ^3.5 xlO 19 W Hz -1 with a linear reso- 
lution of about a 100 pc (Stevens et al. 1999). The increase 
in luminosity for the more distant systems is largely due 
to the coarser linear resolution of the observations. To ex- 
amine this we have convolved a multi-configuration VLA 
image of M82 with a linear resolution of ^6 pc to a reso- 
lution comparable with our observations of NGC6951. In 
the convolved image, the individual components are indis- 
tinguishable and the peak brightness has increased by a 
factor of ^16, corresponding to a luminosity of 6.1 xlO 20 
W Hz -1 at 5 GHz. This is more luminous than the in- 
dividual radio components in NGC6951. A more detailed 
comparison would require observations of NGC6951 with 
higher linear resolution. 

We estimate the SN rate in NGC6951 using the model 
of Condon & Yin (1990) which was based on data for our 
Galaxy. Assuming that the entire emission at 1465 MHz 
(Hummel et al. 1985) is non-thermal, we estimate the SN 
rate to be ~0.07 yr _1 . Bearing in mind that the emission 
at 1465 MHz could also contain a thermal contribution, 
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the SN rate in NGC6951 appears to be reasonably consis- 
tent with estimates of ~0.1 yr^ 1 for M82 (e.g. Huang et 
al. 1994), <0.1 to 0.3 for NGC253 (Ulvestad & Antonucci 
1997), ~0.1 yr _1 for the clumpy irregular starburst galaxy 
Mrk 325 (Condon & Yin 1990) and also for the starburst 
galaxy NGC3448 of the Arp 205 system (Noreau & Kro- 
nberg 1987). 

3.3. The central component 

By fitting a two-dimensional gaussian and an underlying 
tipped plane to the central component, we find it to be re- 
solved with a deconvolved angular size of 0.7x0.2 arcsec 2 
along a PA of 156°. This corresponds to a linear dimension 
of ^80x20 pc 2 . HU also find evidence of a similar exten- 
sion but are unable to confirm it because of the coarser res- 
olution of their observations with the restoring beam being 
elongated along a similar PA. The Ha image of Knapen 
(1996) also seems to show a similar extension. The spectral 
index of the nuclear component obtained by smoothing the 
A3. 6 cm image to that of the lower-resolution A6 cm image 
of Saikia ct al. (1994) is ~0.6. The scaled-array observa- 
tions of HU yield slightly higher values of 0.69 and 0.90 
between A20 and 6 cm for the peak and integrated flux 
densities respectively. The supernova rate estimated from 
the nuclear radio flux density (cf. Colina & Perez-Olea 
1992) is about 0.003 yr" 1 , which is much smaller than the 
SN rate in bright nuclear starbursts. As pointed out by 
Perez et al. (2000), the low SN rate is consistent with the 
faintness of the nucleus in Ha. The resolved central com- 
ponent, which is elongated and has a steep spectral index, 
could be due to a small-scale, weak nuclear jet from the 
Scyfcrt /LINER nucleus. 

4. Conclusions 

High-resolution radio observations of the circumnuclear 
region of NGC6951 reveal a ring of radio emission with 
many discrete components and a slightly resolved nuclear 
component. The discrete components have a median lu- 
minosity of 10 19 W Hz -1 when observed with a linear 
resolution of ~90 pc, and are likely to consist of a mix- 
ture of both SNRs and H n regions. The SN rate has been 
estimated to be ^0.07 yr^ 1 , which is reasonably consis- 
tent with those found in some of the archetypal starburst 
galaxies. The central component, which is resolved and 
has a steep radio spectrum, could be due to a small-scale 
jet from the active nucleus. 
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